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XT] ■ ABSTRACT 

G\ 1 An ever-increasing body of evidence suggests that weakly interacting massive particles (WIMPs) 

constitute the bulk of the matter in the Universe. We illustrate how experimental data, dimensional 
analysis and Standard Model particle physics are sufficient to evaluate and compare the potential of 
detectors searching for such particles either directly (e.g. by their scattering in germanium detectors), 
. or indirectly (e.g. by observing their annihilation into neutrinos in underground detectors). 
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1. Introduction and Results 
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It has become widely accepted that most of our Universe is made of cold dark 
C\| \ matter particles. Big bang cosmology implies that these particles have interactions 

of order the weak scale, i.e. they are WIMPs.S In the early Universe WIMPs are 
in equilibrium with photons. When the Universe cools to temperatures well below 
the mass m x of the WIMP their density is Bolt zmann- suppressed as exp(— m x /T) 
and would, today, be exponentially small if it were not for the expansion of the 
Universe. At some point, as a result of this expansion, WIMPs drop out of equilibrium 
with other particles and a relic abundance persists. The mechanism is analogous to 
nucleosynthesis where the density of helium and other elements is determined by 
competition between the rate of nuclear reactions and the expansion of the Universe. 

For WIMPs to make up a large fraction of the Universe today, i.e. a large fraction 
of Q, their annihilation cross section has to be "just right". The annihilation cross 
^ , section can be dimensionally written as a 2 / m x , where a is the fine-structure constant. 

It then follows that 
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The critical point is that for Q ~ 1 we find that m x ~ m\y, the mass of the weak 
intermediate boson. There is a deep connection between critical cosmological density 
and the weak scale. Weakly interacting particles which constitute the bulk of the 
mass of the Universe remain to be discovered. When our galaxy was formed the cold 
dark matter inevitably clustered with the luminous matter to form a sizeable fraction 
of the 

p x = 0.4 GeV/cm 3 (2) 

galactic matter density implied by observed rotation curves. Unlike the baryons, the 
dissipationless WIMPs fill the galactic halo which is believed to be an isothermal 
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sphere of WIMPs with average velocity 



v x = 300 km/sec . (3) 

In summary, we know everything about these particles (except whether they really 
exist!). We know that their mass is of order of the weak boson mass with: 

tens of GeV < m x < several TeV . (4) 

Lower masses are excluded by accelerator and (in) direct searches while masses beyond 
several TeV are excluded by cosmological considerations. We know that WIMPs 
interact weakly. We also know their density and average velocity in our Galaxy given 
the assumption that they constitute the dominant component of the density of our 
galactic halo as measured by rotation curves. 

Two general techniques, referred to as direct (D) and indirect (ID), are pursued 
to demonstrate the existence of WIMPsB In direct detectors one observes the energy 
deposited when WIMPs elastically scatter off nuclei. The indirect method infers the 
existence of WIMPs from observation of their annihilation products. WIMPs will 
annihilate into neutrinos; massive WIMPs will annihilate into high-energy neutrinos 
which can be detected in high-energy neutrino telescopes. Throughout this paper 
we will assume that such neutrinos are detected in a generic Cherenkov detector 
which measures the direction and, to some extent, the energy of a secondary muon 
produced by a neutrino of WIMP origin in or near the instrument.! It can also detect 
the showers initiated by electron-neutrinos. 

The indirect detection is greatly facilitated by the fact that the sun represents 
a dense and nearby source of accumulated cold dark matter particles.! Galactic 
WIMPs, scattering off nuclei in the sun, lose energy. They may fall below escape 
velocity and be gravitationally trapped. Trapped WIMPs eventually come to equi- 
librium temperature and accumulate near the center of the sun. While the WIMP 
density builds up, their annihilation rate into lighter particles increases until equi- 
librium is achieved where the annihilation rate equals half of the capture rate. The 
sun has thus become a reservoir of WIMPs which we expect to annihilate mostly into 
heavy quarks and, for the heavier WIMPs, into weak bosons. The leptonic decays of 
the heavy quark and weak boson annihilation products turn the sun into a source of 
high-energy neutrinos with energies in the GeV to TeV range. 

The performance of future detectors is determined by the rate of elastic scattering 
of WIMPs in a low-background, germanium detector and, for the indirect method, by 
the flux of solar neutrinos of WIMP origin. Both are a function of WIMP mass and 
of their elastic cross section on nucleons. In standard cosmology WIMP capture and 
annihilation interactions are weak, and we will suggest that, given this constraint, 
dimensional analysis is sufficient to compute the scattering rates in germanium de- 
tectors as well as the neutrino flux from the measured WIMP density in our galactic 
halo. We will derive and compare rates for direct and indirect detection of weakly 
interacting particles with mass m x ~ mw assuming 
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1. that WIMPs represent the major fraction of the measured halo density. Their 
flux is 

0.4 GeV fi cm 1.2 x 10 7 2 , 
d, = n Vy = - 3 x 10° — = cm z s , (5) 

m x cm 6 s ^ X GeV 

where m x c e v = GeV) is in GeV units. 

2. a WIMP-nucleon interaction cross section based on dimensional analysis 

,2 1 



a 



( X N) = (G F m 2 N ) 2 = a DA = 6 x 10" 42 cm 2 , (6) 



3. that WIMPs annihilate 10% of the time in neutrinos (this is just the leptonic 
branching ratio of the final state particles in the dominant annihilation channels 
XX W + W~ or QQ, where Q is a heavy quark). 

Clearly the cross section for the interaction of WIMPs with matter is uncertain. 
Arguments can be invoked to raise or decrease it. Important points are that i) our 
choice represents a typical intermediate value, ii) all our results for event rates scale 
linearly in the cross section and can be easily reinterpreted, and iii) the comparison 
of direct and indirect event rates is independent of the choice. 

We present a simple and totally transparent analysis in which the event rates of 
detectors are derived from the above assumptions. It finesses all detailed dynamics 
and gives answers that are sufficiently accurate considering that the mass of the 
particle has not been pinned down. We will find that the event rate in a direct 
detector is proportional to the WIMP cross section and flux and the density of targets 
m^ 1 , i.e. 

= — <P x *daN(A d ) = (kg)- 1 (year)- 1 , (7) 

dm vtx jv m X GeV 

where represents the number of direct events per unit of target mass. N(A D ) 
represents the coherent enhancement factor for a nuclear target of atomic number 
Ap, e.g. 76 for Germanium, 



N(A) = A s 



I _|_ m x 

m N 



The rates for indirect detection are 

dN m /dA ~ {1.8 x 10- 2 m xGeV } {p(Ajd)N(A id )} {l + 1.9 X 10- 4 m X Gev}" 7 , (9) 

where dNm/dA, in units of (10 4 m 2 ) _1 (year) _1 , represents the number of events from 
the sun per unit area A detected by a neutrino telescope. The factor {pN} should 
be summed over all elements in the sun. Because of additional nuclear form factor 
effects which are neglected in Eq. [] it is adequate to consider oxygen with a solar 
abundance of p = 1.1 % and Aid = 16 as a "typical" element. The observed average 
muon energy should be in the range 1/4 ~ l/6m X G e v- 
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The above parametrizations readily lead to the conclusion that the direct method 
is superior if the WIMP interacts coherently on nuclei (which has been assumed 
for Eqs. and, if its mass is lower or comparable to the weak boson mass my/. 
We will show that in all other cases, i.e. for relatively heavy WIMPs and for all 
WIMPs interacting incoherently, the indirect method is competitive or superior, but 
it is, of course, held hostage to the successful deployment of high energy neutrino 
telescopes with effective area in the ~ 10 4 -10 6 m 2 range and with appropriately low 
threshold. Especially for heavier WIMPs the indirect technique is powerful because 
underground high energy neutrino detectors have been optimized to be sensitive in 
the energy region where the neutrino interaction cross section and the range of the 
muon are large. A kilometer-size detector probes WIMP masses up to the TeV-range, 
beyond which they are excluded by cosmological considerations. 

For high energy neutrinos the muon and neutrino are aligned, with good angular 
resolution, along a direction pointing back to the sun. The number of background 
events of atmospheric neutrino origin in the pixel containing the signal will be small. 
The angular spread of secondary muons from neutrinos coming from the direction of 
the sun is well described by the relation! ~ 1.2°/ ^E^TeV). Measurement of muon 
energy, which may be only up to order of magnitude accuracy in some experiments, 
can be used to infer the WIMP mass from the angular spread of the signal. The 
spread contains information on the neutrino energy and, therefore, the WIMP mass. 
More realistically, measurement of the muon energy can be used to reduce the search 
window around the sun, resulting in a reduced background. 

Before proceeding, we comment on our ansatz for the elastic WIMP-nucleon scat- 
tering cross section. The simplest dimensional analysis implies that the cross section 
is G 2 F m 2 N . This correctly describes the Z-exchange diagram of Fig la, which is of the 
form 
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(m N + m x ) 2 
X 



(a) 

N 



(b) 1 H 



Q 

Fig. 1. Examples of (a) incoherent and (b) coherent WIMP-nucleon interactions. In (b) the gluon 
is a constituent of the target nucleon and Q is a heavy quark. 

For coherent interactions, which we will emphasize throughout this paper, there 
is an additional suppression factor associated with the exchange of the Higgs particle 
with a mass of order of the weak boson mass; see Fig lb. In the specific diagram 
shown the Higgs interacts with the heavy quarks in the gluon condensate associated 
with the nucleon target. It is of the form 

7 . , g 2 rn 2 N m\ 1 ^ 
F H (m N + m x ) 2 m 2 v ' 
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where gn ~ y/Gp m N describes the condensate. Conservatively, we will use the sup- 
pressed WIMP interaction cross section which is appropriate for coherent scattering. 

2. Derivation of Detection rates 

For the case of direct detection the structure of Eq. [7] is transparent .0 For indirect 
detection the number of solar neutrinos of WIMP origin can be calculated in 5 easy 
steps by determining: 

• the capture cross section in the sun, which is given by the product of the number 
of target nucleons in the sun and the elastic scattering cross section 

a Q = f [1.2 x 10 57 ] a DA . (12) 

This includes a focussing factor / given, as usual, by the ratio of kinetic and 
potential energy of the WIMP near the sun. It enhances the capture rate by a 
factor 10. 

• the WIMP flux from the sun which is given by 

0© = x o- /47rd 2 > (13) 
where d — 1 a.u. = 1.5 x 10 13 cm. 

• the actual neutrino flux, which is obtained after inclusion of the branching ratio. 
From (|),(D and ©,(0) 



0, = KT 1 x = 3X10 5 cm- 2 s- 1 . (14) 

m X GeV 

the probability to detect the neutrino,! which is proportional to 

P = poyR^ with 

p = Avagadro # = 6 x 10 23 

cr„ = neutrino interaction cross section = 0.5 x 10~ 38 22,, (GeV) cm 2 
R/j, = muon range = 500 cm ^(GeV) 
or P = 2x 10- 13 m 2 GeV (15) 

Here we assumed the kinematics of the decay chain 

xx -> w + w- 

with E v = \m x (this would be \m x for Q decay) and = \E V = 

finally, dN m /dA = <f> v P = 1.8 x 10~ 6 m xGeV (year)" 1 (m 2 )" 1 (16) 
where dNis/dA represents the number of events from the sun per unit area (m 2 ) 
detected by a neutrino telescope. 
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We can now summarize our results so far by comparing a 10 4 m 2 neutrino detector, 
an area typical of the instruments now being deployed, with a kilogram of hydrogen: 



dN lD /dA = 1.8 x Kr 2 m xGeV (10 4 m 2 ) -1 (year)- ] 
1.4 



dN D /dM = — (kg)" 1 (year)" 1 

^xGeV 

dN B /dM / I0 4 m 2 \ 7.8 x 10 1 
dN m /dA V kg J m 2 GeV 

Direct detection is superior only in the mass range m x < 10 GeV, but this region is, 
arguably, ruled out by previous searches. Indirect detection is the preferred technique. 
This straightforward conclusion may, however, be invalidated when WIMPs interact 
coherently with nuclei and targets other than hydrogen are considered. We discuss 
this next. 



3. Coherent Nuclear Enhancements 

For WIMPs interacting coherently with nuclei in the detector or in the sun, the 
nuclear dependence of the event rates resides in 

• the target density factor m^ 1 in Eq. [7| which is replaced by (Am^^ 1 , 

• the coherent enhancement factor "A 2 ", 

• the nuclear dependence of the cross sections of Eqs. [T^|n| which is obtained by 
the substitution 



incoherent 



n 2 m % m l n i {Am N fm 2 x 



F {m N + m x ) 2 F (Am N + m 



x) 



coherent 

^2 2 



m 2 N ml 1 ^ (9h\ 2 (Am N ) 2 m 2 x 

(m N + m x ) 2 mly \m w J (Am N + m x ) 2 ~ 



The coherent enhancement factor for a nucleus A, including a factor A 1 for 
the target density, is therefore given by 



2 

I _|_ HlX. ~ 

m N 

m N 



1 A 2 (Am N ) 2 m x (m N + m x ) 2 _ ^ (m N + m x ) 2 _ ^ 
A(AmN + m x ) 2 m 2 N m x (ArriN + m x ) 2 

This yields Eq. |[ 

4. Event Rates for WIMPs with Coherent Interactions 

Our simple evaluations, made so far, overestimate the indirect rates for very heavy 
WIMPS because high energy neutrinos, created by annihilation near the core, may 
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be absorbed in the sun. Absorption is stronger for neutrinos and, therefore, mostly 
antineutrinos form the signature for very heavy WIMPS. The probability that an 
antineutrino escapes without absorption is well parametrized by (1 + 3.8 x 10~ 4 E U )~ 7 , 
where E v ~ m x /2. This yields our final result of Eq. |]. 

The relative merits of the two methods are illustrated in the following table, which 
is obtained from Eqs. [^|and establishes that a kilogram of germanium and a 10 4 m 2 
are competitive. 



Table 1. Event rates and signal to background (N/ B). 



m x (GeV) 


Direct (/kg/year) 


Indirect (/10 4 


m 2 /year) 




events Nj B 


events 


N/B 


50 


2.2 x 10 3 7 


2.3 x 10 1 


~ 1 


500 


1.1 x 10 3 7 


2 x 10 2 


~ 10 2 


2000 


2.9 x 10 2 1 


1.7 x 10 2 


~ 10 4 



At the lower energy the event rates for the indirect method are underestimated 
because also the Earth is a source of neutrinos of WIMP origin. 

We conclude that the direct method yields more events for the lower masses, even 
when compared to a 10 6 m 2 detector. As expected, the indirect method is competitive 
for heavier WIMPs with a detection rate growing like El or m 2 . A 10 5 m 2 instrument 
covers the full WIMP mass range, even if the WIMPs do not coherently interact with 
nuclei in the sun. These conclusions are reinforced after considering the signal-to-noise 
for both techniques which we discuss next. 



5. Backgrounds 

Indirect Background. For the indirect detection the background event rate is deter- 
mined by the flux of atmospheric neutrinos in the detector coming from a pixel around 
the sun.B The number of events in a 10 4 m 2 detector is ~ 10 2 /£' A1 (TeV) and the pixel 

size is determined by the angle between muon and neutrino 

m x /4 we obtain 



1.2°, 



£ M (TeV). Using 



the kinematics E^ ~ 



B 



10 2 /^(TeV) 1.1 x 10 5 



ID 



2tt, 



1 2° £ .. 
_ 180 s 



1 2 



m x GeV 



per 10 4 m 2 per year 



This is only valid for large m x , i.e. for E^ 
imation we obtain 



m x /4 > 100 GeV. Without this approx- 



E^GeV) 


# bkgd. events 
in 10 4 m 2 
in 2tt 


# pixels of solar 
size in 2tt 


bkgd. events 
per 10 4 m 2 
per pixel, per year 


10 


3200 


140 


23 


100 


1060 


1.4 x 10 3 


0.8 


1000 


110 


1.4 x 10 4 


8 x 10~ 3 
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For large m x the signal to background ratio is 



N\ dN m /dA » o 

1 - 11 -7.2xl0- 6 m? 



Clearly, the extremely optimistic predictions for signal-to-noise are unlikely to survive 
the realities of experimental physics. One expects, typically, to measure muon energy 
only to order-of-magnitude accuracy in the initial experiments. The energy of showers 
initiated by electron neutrinos should be determined to a factor 2. It is not excluded 
that future, dedicated experiments may do better. The conclusion that high energy 
muons pointing at the sun represents a superb signature, is unlikely to be invalidated. 
Direct Background: about 300 events per year per kg.! Signal-to-noise therefore 
exceeds unity up to 2 TeV WIMP mass. 

These considerations were used to estimate the signal-to-noise N/B in Table 1. 



6. Dynamics? 



We emphasize that our results are representative for the specific and much stud- 
ied example where the lightest supersymmetric particle is Nature's WIMP.Q Clearly 
dynamics, which is now defined, can alter our conclusions, but only in "conspirato- 
rial" ways. Dynamics can, on the other hand, increase rates as well, sometimes by 
well over an order of magnitude, over and above the rates obtained from dimensional 
analysis in this paper. Our qualitative conclusions are valid, at least in some average 
sense, in supersymmetry. Our results do, in fact, closely trace the supersymmetry 

prediction of reference 2 for the choice of Higgs coupling = 1, in their notation. 

We feel that the development of detectors should be guided by an analysis like 
ours rather than by dynamics of theories beyond the standard model for which there 
is, at present, no experimental guidance. 
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